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On the penultimate tail behavior of Weibull-type models 

Abstract: The Gumbel max-domain of attraction corresponds to a null tail index which do not distin- 
guish the different tail weights that might exist between distributions within this class. The Weibull-type 
distributions form an important subgroup of this latter and includes the so-called Weibull-tail coefficient, 
usually denoted 0, that specifies the tail behavior, with larger values indicating slower tail decay. Here we 
shall see that the Weibull-type distributions present a penultimate tail behavior Frechet if 6 > 1 and a 
penultimate tail behavior Weibull whenever 6 < 1. 
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1 Introduction 

The main objective of an extreme value analysis is to estimate the probability of events that are 
more extreme than any that have already been observed. By way of example, suppose that a 
sea-wall projection requires a coastal defense from all sea-levels, for the next 100 years. The use of 
extremal models enables extrapolations of this type. The central result in classical Extreme Value 
Theory (EVT) states that, for an i.i.d. sequence, {X„} n >i, having common distribution function 
(d.f.) F, if there are real constants a n > and b n such that, 

P(maxpTi, ...,X n ) < a n x + b n ) = F n (a n x + b n ) — > n ^oo G y (x) , (1) 

for some non degenerate function G 7 , then it must be the Generalized Extreme Value function 
(GEV), 

G-y{x) =exp(-(l + 7.t)- 1/7 ), 1 +jx > 0, 7 G K, 

(Gq(x) = exp(— e~ x )) and we say that F belongs to the max-domain of attraction of G 7 , in short, 
F G T)(G~f). The parameter 7, known as the tail index, is a shape parameter as it determines the 
tail behavior of F, being so a crucial issue in EVT. More precisely, if 7 > we are in the domain of 
attraction Frechet corresponding to a heavy tail, 7 < indicates the Weibull domain of attraction 
of light tails and 7 = means a Gumbel domain of attraction and an exponential tail. 

However, as Fisher and Tippett ([1], 1928) remarked, if one approximates the distribution of 
the successive maxima of normal samples not by the limit distribution Gumbel but by a sequence 
of other extreme value distributions converging to the limit distribution, the approximation is 
asymptotically improved. They called penultimate distributions to this sequence of approximating 
extreme value distributions. 

Here we will analyze the penultimate tail behavior of the Weibull-type models which are a wide 
class in the Gumbel max-domain. 

Weibull-type models have representation 

1 - F(x) = exp(-H(x)), x>x o >0,e>0 (2) 

where 

H(x) = x 1 ' 6 ^) or H~ l {x) = x°l*(x), (3) 

with H^ 1 denoting the generalized inverse of H and functions I and I* are slowly varying at infinity 
(i.e., l(tx)/l(t) — > 1 as t — > 00 for all x > and the same holds for I*). 

We say that functions H and are regularly varying with indexes, 1/0 and 6, respectively. 

The parameter 9, called the Weibull-tail coefficient, governs the tail behavior of F, with larger 
values indicating slower tail decay. The Weibull-type distributions form an important subgroup 
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within the Gumbel class (7 = 0) an the tail behavior can then be specified using the Weibull-tail 
coefficient (Dierckx et ai, [5] 2009). 

Weibull-type models include well-known distributions such as, Normal (9 — 1/2), Weibull(a,A) 
(9 = I /a), Extended Weibull(/3,<5) (9 = 1//3), Exponential, Gamma and Logistic (9 = 1) (Gardes 
and Girard [5J, 2008). 

We will prove that the Weibull-type models present penultimate tail behavior Frechet or Weibull 
whenever > 1 or 9 < 1, respectively. To this end, we will use a result in de Haan and Gomes 
([5], 1999) and in Gomes (0, 1984). 

We remark that the case 9=1 cannot be treated in a unified way since we can find penultimate 
tail indexes with different orders. For instance, the Exponential and Logistic have j n of order l/n 
but the Gamma distribution have j n of order 1/ log 2 n (Gomes [7|, 1993). Therefore, we shall 
assume that 9 ^ 1 all over the paper. 

2 Results 

Consider 

k{x) := ±[-\og{-\ogF{x))} (4) 



Observe that 
and, by ©, 
leading to 

Therefore 



-\ogF{x) = (1- F(x))(l+0(1- F(x))) 
log(- \ogF{x)) = H(x) - log(l + 0(e- H ^)), 

0(e- H ^ 



k(x) = H'{x)(l + 1 + g (e _ g( i )} ) = H'(x)(l + o(l)). (5) 



(0) 



w= *"(*)(* + i^gay) + p^)] 2 0(e ' 2flW) aS;:j!;; o(e ' flW)) 

= H"{x){l+o{\)). 
Analogously we derive 

k"{x) = H"'{x)(l + o(l)) and k'"{x) = H {m \x){\ + o(l)). (7) 
Consider H(x) given in ([3]). We have 

H'(x)= X ^-H( X )^ + ^) 

H"(x) = x^-Hix)^- 1 ^- 1 - 1) + 20-^ + 2^2) 

= x rl - 3 l{x)(9- 1 {9- 1 - l)^- 1 - 2) + Stf- 1 ^ 1 - l)2g£l + ZO' 1 ^^- + ^gp) (8) 

i/(™)(.T) = a; 8 " 1 -^^)^- 1 ^- 1 - l)(9- 1 - 2){e~ 1 - 3) + 46* 1 (6* 1 - l)(9- 1 - 2)2^1 
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Assuming that, as x — > oo, 

T^-^O. ^^0, £!g£)^Oand 2^ 0, (9) 
and applying (HJ-Q, we have that 



xk(x) = 6>-V _1 Z(a:)(l + o(l)) 

x 2 /c'(a;) = xk^id- 1 - 1) 

x 3 /c"(a;) = xkix^e- 1 - 2)(6» 1 - 1) 

x A k"'{x) = xkix)^- 1 - 3){0- 1 - 2){9- 1 - 1) 



(10) 



Note that, if l(x) is monotone for x > x$ > 0, then xl'(x)/l(x) — > 0, as x — > oo. The other 
conditions in © are also satisfied by the most common models. 

Observe that equations ©-([7]) hold for the class of d.f.'s 

- log F{x) = cxp(-H{x)) (11) 

which have been studied in Canto e Castro ([2], 1992). More precisely, it was derived the structure 
of the remainder F n (a n x + b n ) — G 1 (x), by showing that F is in class Ai of Anderson ([1], 1976), 
i.e., 

Um k "^ = Um x3k "^ = Um i^ 1 - 2 ) = o (12) 

x^-oc k(x)k'(x) x^ao xk{x)x 2 k'{x) a->oo xk(x) 



leading to 



F n (a n x + b n ) - G,(x) = y (^l)p 7 (x)(l + (1)), (13) 
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uniformly for a; in bounded intervals in the support of G 7 , where g 7 (x) = G' y (x), a n — l/k(b n ) 
and F{b n ) = exp(— 1/n). More details can be seen in Canto e Castro ([2J, 1992). 

Observe that 

b n = H- 1 (logn), (14) 
and the rate of convergence of F n (a n x + b n ) to G 1 (x) is 

k'(b n ) b 2 n k'{b n ) 6- 1 -! 1-6 



k 2 (b n ) {b n k(b n )) 2 fl-iflogn) logn 



(15) 



Now consider x F := supjx : F(x) < 1} and 

ip{t) = (l/fc)'(t) = -fc'(f)/(fc(t)) 2 . (16) 
Based on the von Mises' first order condition, 

lim p(t) = 0, (17) 

t—>x F 

the von Mises' second order condition, 

lim , f.® . = 0, (18) 
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and von Mises' type penultimate condition, 

lim , f.® . = 0, (19) 

Theorem 1 in Haan and Gomes ([8], 1999) allow us to derive bounds for (F n (a n x+b n )—G 7n (x)) / j' (log n) , 
with a n and b n given above, where 

M-^ 1 ™ <20) 

and the penultimate tail index is 

7n = 7 (logn) - -k'(b n )/k 2 {b n ). (21) 

Proposition 2.1. Weibull-type models satisfy the von Mises' conditions il7\) - HT^) whenever con- 
ditions in fPP are fulfilled. Moreover, they present penultimate tail behavior Frechet if 9 > 1 and 
Weibull if9<l. 

Dem. First observe that x F = +oo, 

k"(t)k(t)~2(k'(t)f 
Wt) = WW 

and 

Now we have, successively, 



v 1 ' (k(tW V k(t) V kit) ) ) (k(tW 



, n ,. k'(t) ,. t 2 k'(t) ,. ^-l l + ol 

lim = lim - A „ = 1 m - - ,\' = lim — \ a . V^r = 0, (22) 

t^^ K ' t-+x* (k(t)) 2 t^x*- (tk(t)) 2 t^xF 6-H0- 1 (t)(l + o 1)) ' 1 ; 



by (d2) and (J22J), 



and, applying ((T 



V (*) ,. A, * (*) fc (*) \ „ 

lim tttHt = lim (2 ., /.(„ - , , .;/, J = 0, 23) 
t^x*- k(t)ip{t) 2 k(t)k'(t)J ' v ; 



2k'{H- 1 {t)) 2 - fc ,, (ff- 1 (t))fc(g- 1 ft)) 

and applying (jT4j we have 



(24) 



6k'(t)k"(t)k(t) - 6(k'(t)) 3 - k"'(t)(k(t)) 2 
J£* k(t)<p'(t) ~ A (k(t)) 2 (2(k'(t)) 2 - k"(t)k(t)) 

6(6>- 1 - l)^- 1 - 2) - 6(6>- 1 - l) 2 - (g- 1 - 2)(6'- 1 - 3) 
t-J£* t/fc(t)[2(0- 1 -l)-(0- 1 -2)] 

By (|15p . the penultimate tail index, 7„, in (|2ip is given by 

ln = 7 (Iogn) = -k'(b n )/k 2 (b n ) ~ ^i. (25) 

Therefore, we obtain 7„ > and 7„ < if, respectively, > 1 and 6> < 1. 

Now we compute the rate of convergence 7'(logn). Observe that, after some calculations, 



7 (*) = ^rr-l W u ' ( 26 ) 



x 2(b n k(b n )(6-^ - l)) 2 - (b n k(b n )) 2 (9-i - - 2) 20(1-0) 

7(lognH (MKM? lbg^F D (27) 



Remark 2.1. The result in Provosition \2.1\ can also be derived based on Gomes (f^jj, 1984)- This 
latter imposes the condition 

lim , . f. ^,,„ = c < oo (28) 

where ip(t) = (l/k)'(t) = —k'{t)/(k(t)) 2 , and obtains boundedness of (F n (a n x + b n ) — G 7n (x))/"f 2 , 
with penultimate tail index, "f n , given in 1125]) . 
Observe that 

<p>(t) _ 2 k"(t)k(t) 



k{t){^{t)Y (k'(t)) 2 ■ 

Applying UP]) , we obtain 

k"{t)k(t) _ t 3 k"(t)tk(t) _ (Q- 2 - W- 1 +2 + o(l))(l + o(l)) _ (6- 1 - 2){6- 1 - 1) 
(k'(t)) 2 (t 2 k'(t)) 2 ~t^%* (9- 1 -l + o(l)) 2 ~ (6*- 1 - l) 2 

and hence 

lim ™ - 1 



t->x^ k(t)(ip(t)) 2 1-6' 

Thus condition 128]) holds for 8 =/= 1 (the situation that we are always considering; see the last 
paragraph of section^. The rate of convergence is given by j 2 = (7(logrt)) 2 = (k' (b n ) / k 2 (b n )) 2 = 
((8 — l)/logn) 2 , i.e., also the same order obtained in HFfy . 
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